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We calculate the back reaction of cosmological perturbations on a general relativistic variable which 
measures the local expansion rate of the Universe. Specifically, we consider a cosmological model in 
which matter is described by a single field. We analyze back reaction both in a matter dominated 
Universe and in a phase of scalar field-driven chaotic inflation. In both cases, we find that the 
leading infrared terms contributing to the back reaction vanish when the local expansion rate is 
measured at a fixed value of the matter field which is used as a clock, whereas they do not appear 
to vanish if the expansion rate is evaluated at a fixed value of the background time. We discuss 
possible implications for more realistic models with a more complicated matter sector. 



I. INTRODUCTION 

Because of the nonlinear nature of the Einstein equa- 
tions, linear fluctuations about a homogeneous and 
isotropic FRW cosmology will back-react on the back- 
ground on which they live. The back-reaction of short 
wavelength gravitational waves was studied a long time 
ago by Brill, Hartle and Isaacson Q). In this approach, 
the back-reaction of fluctuations on the spatially aver- 
aged metric can be described by an effective energy- 
momentum tensor which contains the spatially averaged 
terms in the Einstein tensor which are second order in 
the amplitude of the fluctuations. More recently, Tsamis 
and Woodard Q have initiated a detailed study of the 
back-reaction of long wavelength gravitational waves in 
a de Sitter background. Scalar-type metric fluctuations 
(called "cosmological fluctuations" in what follows) are 
also expected to back-react on the background space- 
time. Since in the context of inflationary and post- 
inflationary cosmology the scalar metric fluctuations are 
believed to dominate over the effects of gravitational 
waves, the back-reaction of these cosmological perturba- 
tions are expected to give the dominating back-reaction 
effect. 

In the formalism of Q was generalized to de- 

scribe the back-reaction of cosmological perturbations on 
the spatially averaged metric. On this basis, it was ar- 
gued [^|-^| that gravitational back-reaction in scalar field- 
driven inflationary models, calculated up to quadratic 
order in perturbations and to leading order in the long 
wavelength expansion and in the slow roll approximation 
could decrease the expansion rate of the universe and po- 
tentially solve the cosmological constant problem (see 
also for similar results in the context of back-reaction 



studies of gravitational waves). The formalism of 
covariant under first order space-time diffeomorphisms []. 

However, as emphasized by Unruh [M [j], the approach 
of H||,H is deficient in several respects. First of all, due 
to the nonlinear nature of the Einstein equations, calcu- 
lating an "observable" from the spatially averaged metric 
will not in general give the same result as calculating the 
spatially averaged value of the observable. More impor- 
tantly, the spatially averaged metric is not a local physi- 
cal observable. Thus, to take into account the deficiencies 
of the previous work on gravitational back-reaction, we 
must identify a local physical variable which describes 
the expansion rate of the Universe, calculate the back- 
reaction of cosmological perturbations on this quantity, 
and then take the spatially averaged value of the result. 
It is important to fix the hypersurface of averaging by a 
clear physical prescription in order to remove the possi- 
bility of being misled by effects which are second order 
gauge artifacts. 

In this paper, we propose an implementation of this ap- 
proach. We focus on the variable which yields the general 



*The back-reaction of small-scale (i.e. smaller than the Hub- 
ble radius) cosmological perturbations has been considered 
(without taking into account the issues of gauge freedom) 
in p[-p^[. This problem has also been considered in the con- 
text of Newtonian cosmological perturbation theory in 
More recently, Nambu ]l8|-pc| has initiated a program to com- 
pute back-reaction effects on the spatially averaged metric 
using the renormalization group method. 

'We are also grateful to Andrei Linde and Alan Guth for 
detailed private discussions on these points. 
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relativistic definition of the local expansion rate, calcu- 
late this quantity to second order in the amplitude of 
the cosmological fluctuations, in terms of a time variable 
defined by an unambiguous physical prescription. For 
simplicity, we assume that matter is described by a sin- 
gle field (either a hydrodynamical field or a single scalar 
field). 

We study two examples, first a matter dominated Uni- 
verse, and second the inflationary phase of a cosmol- 
ogy dominated by a single scalar field. In both cases 
we find that the leading infrared contributions to the 
back-reaction on the local expansion rate of the Universe 
vanish, in contrast to the findings of the initial work on 
gravitational back-reaction of cosmological fluctuations 
and confirming the analysis of We thus con- 
firm the conclusions reached recently in p2p3| ] where a 
different variable related to the local expansion rate is 
proposed, and different techniques to evaluate this vari- 
able are used. 

Note that when evaluated at a fixed background time, 
our leading infrared back-reaction terms give a non- 
vanishing contribution. This leads us to the conjecture 
that in more realistic models in which a second field is 
present to determine time (e.g. the microwave back- 
ground), the leading infrared back-reaction terms will not 
vanish. 



II. A LOCAL OBSERVABLE 

For a general perfect fluid flow in a curved space-time 
we consider the velocity four-vector field u a tangential 
to a family of world lines. In the context of cosmology, 
we can always define a preferred family of world lines 
representing the motion of a set of comoving observers. 
In the case of hydrodynamical matter, this is easy since 
the energy-momentum tensor is already defined in terms 
of a velocity four-vector field. Also in the case of scalar 
field matter, a corresponding four-vector field can be de- 
fined, although a bit more care is required to obtain a 
consistent definition. In both cases we have 



is the local expansion rate of the tangential surfaces 
orthogonal to the fluid flow, u> is the vorticity tensor 
(with uj a pvr — 0), and a is the shear tensor (satisfying 
o a pvP = 0, u a a = 0). 

For a homogeneous Universe with scale factor a(t) the 
Hubble expansion rate H is 



H = a/a = -6 . 
' 3 



(5) 



For a cosmological model with fluctuations, is local 
in space and time. We will use to define the local 
expansion rate S/S via the local analog of (^), namely 
via the equation 



S/S 



(6) 



The quantity S/S is a much better measure of the locally 
measured expansion rate in a Universe with fluctuations 
than the Hubble expansion rate used in HQ determined 
via the spatially averaged metric including back-reaction. 
It is a mathematically simpler object than the variable 
recently introduced in p2| which involves the integral 
along the past light cone of the observation point. If 
we are interested in evaluating the expansion rate for a 
typical observer, we propose to take the spatial average 
of the local expansion rate defined via @ . 

Now that we have defined the observable we are inter- 
ested in, the procedure will be as follows. First, we must 
determine the velocity four-vector field u a for the models 
we are interested in. Then, we use the Einstein equations 
to express u a in terms of the metric perturbation. Tak- 
ing the relative amplitude of the metric fluctuations as 
the expansion parameter, we then calculate 0, our local 
measure of the Hubble expansion rate, to second order. 
After evaluating the result on a physically determined 
hypersurface we can then study the back-reaction of cos- 
mological fluctuations on the locally measured Hubble 
expansion rate. In this paper we will focus on the leading 
infrared contributions to back-reaction, the terms found 
to dominate the back-reaction effects in 



u a u a = 1 



(1) 



where a runs over the space-time indices. The projection 
tensor onto tangential three-surfaces orthogonal to u a is 



h a (3 = 9a[3 — U a Up . 



(2) 



The first covariant derivative of the four-velocity can 
be decomposed as (see e.g. |24] for details ) 

U a ;f3 = + Cr aj3 + -Qh a j3 - U a 1lp . (3) 



Here 



(4) 



III. DERIVING THE EXPANSION RATE FOR 
SCALAR METRIC PERTURBATIONS 

In this section we consider a model with hydrodynam- 
ical matter. Starting from the expression for the metric 
to linear order in the fluctuations <E> (see |^5| for a de- 
tailed review), we determine the velocity four vector field 
u a to the order required to analyze the leading infrared 
terms in the back-reaction to quadratic order. To obtain 
the full back-reaction terms (including terms which dom- 
inate in the ultraviolet but are negligible in the infrared) 
we should calculate u a consistently up to second order. 
However, if we are only interested in the leading infrared 
terms, it is sufficient to keep all the terms quadratic in $ 
but not containing any spatial gradients. 
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In order to obtain the complete result for gravitational 
back-reaction we would have to look at the Einstein equa- 
tions for a perfect fluid with energy density p and pres- 
sure P, 



Gftu = (P + p)ujj,u v - Pg^y 



(7) 



(in units in which 8nG = 1), which, since G£J = R, will 
yield 



R = p-3P 



P : 



(8) 
(9) 



and lead to an equation that can be solved perturbatively 
to any desired order for uf. 



(10) 



However, as mentioned above, here we just use the results 
for u l which are of linear order. Since we will calculate 
the divergence of u^, our prescription implies that we are 
ignoring some of the extra second order gradient terms . 
The result can also be used for scalar fields if we define 
the u vector field in a proper way. 

For an unperturbed Robertson- Walker metric, the 
four- velocity field u in comoving coordinates would be 



= (1,0,0,0) 



(11) 



In linear perturbation theory, and in the case of simple 
forms of matter (such as a single fluid or a single scalar 
field) for which there is to linear order no anisotropic 
stress, the metric (in longitudinal gauge) can be written 
as 



ds 2 

1*3 



(?7) 2 ((1 + 2<S>)d v 2 - (1 - 2^) lij dx i dx j ), (12) 
' ' 2 1 - 2 )] (13) 



+ -)C(x 



I) 



where K, = 0,1,-1 depending on whether the three- 
dimensional space corresponding to the hypersurface t — 
const, is flat, closed or open. In this paper we will take 
it to be zero in order to simplify the calculations. The 
time variable 77 appearing in ( |l2| ) is conformal time and 
is related to the physical time t via drj = a~ 1 dt. For the 
forms of matter considered here, 5* = $ at linear order 
0. As discussed e.g. in [p5[ , in longitudinal gauge the 
spatial components of the four-velocity vector field are 
related to $ via 



\H 2 -W + /C)" 1 (a$) / i 



(14) 



"'"Even if we did not make the assumption $ = it turns out 
that at second order all infrared terms depending on ^ — $ 
will drop out. 



where a prime denotes differentiation with respect to r\ 
and TL = a' /a. Using equation ([l]) we can derive the 
expression for the time component of u a in terms of $: 



-<S> 2 



(15) 



Now that we have all components of u a , we take the 
covariant derivative of it and retain all $ dependence up 
to second order ^ These second order terms principally 
come from the Christoffel symbols. Other second order 
terms (which as mentioned before are gradient terms) 
could be added to the ones computed here if we were 
to solve the Einstein equations in the form (|l7J) beyond 
linear order. A straightforward calculation yields 



e = 3— (1-$ + 



-$ 2 ) 

2 ; a 



a 2 (H'-H 2 ) 



(16) 



Since a/a = a 1 / a 2 , we can immediately read off the extra 
terms contributing to the local expansion rate which re- 
sult from the presence of cosmological fluctuations. Upon 
spatial averaging at a fixed conformal time, the terms lin- 
ear in <!> drop out. Hence, it follows that if evaluated at 
a fixed conformal time, infrared modes on average lead 
to an increase in the expansion rate compared to what 
would be obtained at the same conformal time in the ab- 
sence of metric fluctuations. Whether this is a physically 
measurable effect from an observational point of view will 
be discussed in more depth in following sections. 



IV. EXPANSION RATE FOR A MATTER 
DOMINATED UNIVERSE 

Now let us use the result of the last section to (as an 
example) calculate the local Hubble expansion rate for a 
matter-dominated Universe. In this case, the scale factor 
a and the scalar metric perturbation $ have the following 
dependence on the conformal time r\ : 



a(r]) = a m r) 2 jl 
®(ri,x) = d(x) +C 2 (x)r)- 



(17) 
(18) 



where constant, and C\ and C2 are time- 

independent. The second equation is valid in the long 



^This is in the philosophy of the general back-reaction ap- 
proach in which it is assumed that the fluctuations of the 
metric and matter satisfy the linear perturbation equations, 
and we compute their back-reaction on physical quantities to 
second order. It is not a consistent second order perturbative 
formalism. 
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wavelength (super-Hubble-scale) limit and has been ex- 
plicitly derived in |p5| . Now using equations (^6|), ( |l7| ) 
and ( Jlq ) we obtain O in terms of conformal time: 

Q = —r"{i-c 1 {x) + i -ci) 

+ — 7 7 - 6 (i(a i Ci(a;))(a i C2(a;)) + d 2 C 2 (x)) (19) 

+ —^ 8 (3C 2 (x) +36C 1 {x)C 2 {x)) 

+ — V - n (diC2(x)) 2 + — V - 13 (C 2 (x)) 2 

Some of the terms are decreasing very fast as a function 
of time and thus we can ignore them. If we just keep the 
terms with the powers —3 and —1 of r\ then we get: 

e = — »r 3 (i - ci(aO + |c?) (20) 

-\— V-'mC^xjf + dfC^x)} 

If we take the average of on a constant 77 hypersur- 
face, only the terms quadratic in the fluctuation variables 
survive. Thus, considering large values of r\ and focusing 
on the second order terms, it appears from (^0|) that in- 
frared modes give a positive contribution to and thus 
lead to a speeding-up of the expansion, whereas ultravi- 
olet modes enter with a negative sign and thus yield a 
slowing effect, the latter becoming more significant (rela- 
tive to the unperturbed expansion rate) for larger values 
of rj. 

However, before drawing definite physical conclusions 
from our analysis, we must take into account that the 
background time 77 is not an observable quantity. To ob- 
tain results for back-reaction of any real physical signifi- 
cance we have to find an observable variable like proper 
time and evaluate the expansion rate in terms of this 
variable, so that we can discuss its evolution from an 
observer's point of view. 

If we use equation ( |l2| ) for the metric, equation ( |l8| ) for 
<&, and equation (|l7]) for the scale factor, we can find the 
expression for the proper time r in terms of conformal 
time. Since 

dr 2 = a(r/) 2 (l + 2$)d?7 2 , (21) 
a simple integration yields 

r=°f{l + C l -\c 2 W (22) 
-^(C 2 ~C 1 C 2 ) V - 2 + ^C 2 r 1 ^ 

In the approximation of large values of r\ we can ignore 
the second and the third term of this equation and thus 
obtain: 



„-a = ^ (1 + Cl _l C?)r -i. (23 ) 

Now we can use this relation and substitute into the first 
term of equation ( f20|) . We see the effects of the dominant 
infrared terms on the local expansion rate cancel exactly 
up to second order in perturbations, when evaluating 
at a fixed proper time: 

&IR = 2T- 1 . (24) 

This implies that at least in the approximation of keeping 
only the leading infrared terms, there is no local gravita- 
tional back-reaction of cosmological fluctuations on the 
local Hubble expansion rate in this matter-dominated 
universe. 



V. EXPANSION RATE IN TERMS OF SCALAR 
FIELD AS AN OBSERVABLE 

We now move on to the example more relevant to the 
work of namely a Universe dominated by a single 

real scalar field ip, which is a toy model for inflationary 
cosmology. During inflation, fluctuations which are gen- 
erated on sub-Hubble scales early on during the inflation- 
ary phase are red-shifted to scales much larger than the 
Hubble radius. Thus, in this context it is of great interest 
to consider the back-reaction of infrared modes. 

In the following, we will generalize the previous analy- 
sis to be applicable to matter consisting of a single scalar 
field. In this case one can treat the scalar field as a per- 
fect fluid and derive the velocity four-vector field. To do 
this, we need to write the energy-momentum tensor T^ v 
of the scalar field (right hand side of the following equa- 
tion) in the form of an energy-momentum tensor for a 
perfect fluid (left hand side of the following equation): 

{p + P)u^u v - Pg^ v = d^pd^p - C .g^ (25) 

At the level of the background fields, the two expressions 
are identical if we take P = C and u M = Ad^tp with 

A = (d^d u ^y 1/2 . (26) 

Now that we have shown that the energy-momentum 
tensor of a scalar field can be written in the form of that 
of a perfect fluid, we can use the expression (|l^) which 
gives in terms of the metric fluctuations to evaluate the 
local effect of gravitational back-reaction of cosmological 
fluctuations. Let us first for convenience rewrite of 
( |l6| ) in terms of the physical time t: 

e.Ai- $ + V)-3H»!)iM (27) 

a 2 ad — a 1 

Theoretically, the scalar field <p is an observable. In 
fact, in a system with a single matter field <p, it is this 
field which must be used as a clock. Hence, to obtain 
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results with physical meaning, we must evaluate O on a 
surface of constant ip and not constant t. 

As discussed above, in the context of inflationary cos- 
mology it is important to study the effects of infrared 
modes up to second order. We will also assume that we 
are in the slow rolling regime of inflation. In our case 
(as well as for a more general case), the prescription is 
to calculate t, a and in terms of ip, and to insert 
the results into the general expression ( p7| ) for 0. The 
relation between t and ip can be derived starting from: 



e = V3y/v(<p) 



(33) 



ip{t) = <p Q (t) + 5<pi(t) , 
which can be written as: 

t = f^if - Sl Pi(t)} ■ 
Thus, the equation 



(28) 



(29) 



dtp 
2 dip 2 



dt 



(30) 



relates t and (p. 

We now wish to express the value of © in terms of ip 
(note that we are considering the local value of O in this 
analysis, and there is no need to perform a spatial av- 
eraging). In order to relate the metric, its fluctuations 
and the scale factor to the scalar field and its fluctua- 
tions we need to take an explicit form for the potential 
and make use of the Einstein constraint equations. The 
G 0t and G 00 equations relate <£> to <p (at the level of the 
first order fluctuations) and the Hubble parameter of the 
background to the background scalar field ipo (at the level 
of the unperturbed Friedmann equations). In our simple 
system we do not need to go through these calculations 
explicitly (as long as we are interested only in the lead- 
ing infrared terms) since by dropping the gradient terms 
from the G 0i and G 00 equations, it can be shown that 



H 



1 



VI + 2$ V3 



(31) 



If we expand the right hand side in terms of $ we get the 
result 
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h(i - $ + -$ 2 ) = -^vni) 



(32) 



and substituting this result into equation ( J27| ) and ne- 
glecting terms containing $ and di<& (which are sub- 
dominant compared to other terms in the infrared and 
slow roll limit) leads to the final result for the local ex- 
pansion rate: 



which as a function of tp is the same as the relation for 
an unperturbed background. 

Thus, again we do not see any back-reaction of cosmo- 
logical perturbations on the local expansion rate in this 
approximation. In retrospect it is easy to understand this 
result, since by neglecting gradient terms we have a 
Friedmann Universe whose expansion rate satisfies equa- 
tion (p3 



VI. CONCLUSIONS 

In this paper, we have studied the back-reaction effects 
on a local observable which measures the local expan- 
sion rate of the Universe. The observable gives the rate 
at which neighboring comoving observers separate and 
coincides with the usual definition of expansion in the 
context of the fluid approach to cosmology. In order to 
obtain a physical quantity, we evaluated the observable 
at a fixed value of the scalar field. 

We evaluated our observable, the local physical expan- 
sion rate, in a simple toy model of chaotic inflation con- 
sisting of a single scalar matter field coupled to grav- 
ity. We found that the leading infrared terms, the terms 
which dominate the effects discussed in |3] and fa] , cancel 
if we evaluate the observable at a fixed value of the scalar 
field, whereas they do not vanish if we evaluate them at 
a fixed value of the background time. The former result 
is a physical result since it corresponds to a physical ob- 
servable evaluated at a space-time point specified by a 
physical prescription, whereas the latter result (obtained 
by evaluating at a fixed background time) does not have 
a diffeomorphism-invariant meaning. Our analysis thus 
confirms the concern of [|T) that the results obtained in 
H and H are not invariant under second order gauge 
transformations. Our results confirm the conclusions of 
p3f reached by means of a different method of analysis 
applied to a different physical observable. 

Our result does not imply that there is no back- 
reaction of the infrared modes of cosmological pertur- 
bations. There is no reason to expect that the next to 
leading infrared terms in our result will cancel (they do 
not cancel in the analysis of |23) ) . One of the advantages 
of our technique is that they can be evaluated without 
too much trouble. This is left to a future publication. 
So, even in single field models of inflation there might be 
some non- vanishing back-reaction of infrared modes. 

We expect that back-reaction of infrared modes will 
be much more important in two field models of infla- 
tion. Let us assume that the matter sector of the theory 
contains both an inflation field ip and a regular matter 
field x (with nonvanishing and time-dependent spatial 
average) which, for example, could represent the cosmic 
microwave background. In this case, it is no longer true 
that long wavelength fluctuations have no physical effects 
on local observables. If the measurement point is (in an 
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unambiguous physical way) determined by a fixed value 
of the field x, then the local expansion rate may sensi- 
tively depend on the amplitude of the long wavelength 
fluctuations in ip. Thus, the leading infrared terms may 
not cancel when evaluated according to the abovemen- 
tioned physical prescription in the same way that they 
do not cancel in the analysis of this paper when the ob- 
servable is evaluated at a fixed value of the background 
time. 

There is a close analogy with the analysis of the para- 
metric amplification of super-Hubble-scale cosmological 
fluctuations during inflationary reheating. From the 
point of view of the background space-time coordinates, 
it appears pjj ] that the parametric amplification of mat- 
ter fluctuations on super-Hubble scales in an unperturbed 
cosmological background (see e.g. |E8 29 for a discus- 



sion of parametric resonance during reheating) would 
imply the parametric amplification of the cosmological 
fluctuations on these scales. However, it can be shown 
that in single field models physical observables measur- 
ing the amplitude of cosmological fluctuations do not feel 
any resonance 30 32 23]. In contrast, in two field mod- 



els of inflation there is [|33 34 parametric amplification 
of super-Hubble-scale cosmological fluctuations. In this 
case, there is a fluctuation mode corresponding to en- 
tropy fluctuations which cannot locally be gauged away. 
This mode is (in certain theories) parametrically ampli- 
fied during reheating, and in turn drives the parametric 
resonance of the super-Hubble scale curvature fluctua- 
tions. 

Our techniques allow us to calculate the back-reaction 
of cosmological fluctuations in two field models in a very 
similar way to what is presented here. Results will be 
presented in a followup publication |35f| . 
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